New analytic models of "traversable" wormholes 



October 2, 2008 



Alexander Shatskiy^, I.D. Novikov^'^'^, N.S. Kardashev^ 

'66 

§ ! -"^ Astro Space Center, Lebedev Physical Institute, Moscow, Russia, shatskiy@asc.rssi.ru 
CN I ^Niels Bohr Institute, Blegdamsvej 17, DK-2100 Copenhagen, Denmark 
^ • ^Kurchatov Institute, Moscow, Russia 

o 



Abstract 



The analytic solution of the general relativity equations for spherically symmetric 

wo 
all( 



O , wormholes are given. We investigate the special case of a "traversable" wormhole i.e., one 

' allowing the signal to pass through it. The energy- momentum tensor of wormhole matter 



is represented as a superposition of a spherically symmetric magnetic field and dust matter 
with negative matter density. The dynamics of the model are investigated. We discuss 
both the solution of the equation with a A-term and without it. Superposing enough dust 
^ ' matter, a magnetic field, and a A-term can produce a static solution, which turns out to 

. be a spherical Multiverse model with an infinite number of wormholes connected spherical 

j. I universes. Corresponding solution can be static and dynamic. 
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1 Introduction 



A wormhole (WH) [T]-[l9]) is a hypothetical object described by a nonsingular solution of the 
Einstein equations with two large (or infinite) space-time regions connected by a throat. The 
^ I two large space-time regions can be located in one universe or belong to different universes in the 
d I Multiverse model (see [IE])- In the last case, "traversable" WHs afford a unique opportunity 
to explore other universes. 

In the present paper, we analytically study the dynamics of a spherical model of a non- 
equilibrium WH filled with matter. This matter consists of a magnetic field and dust with 
negative mass density. The obtained solution therefore generalizes the Tolman solution [20J for 
a model with a spherically symmetric electromagnetic field. As we see below, this generalization 
leads to essentially new and important solutions. We use the method of calculations of physical 
quantities in the frame co moving with dust. 

As the initial model for a WH, we use a static model in which gravitational accelerations are 
everywhere identically zero. Hence, the effective masses of both WH mouths vanish, although 
the geometry of three-dimensional space is certainly non- Euclidean. Such a model is considered 
in [TQ] and [T7], where all matter is represented by a gravitating scalar field. We change the 
scalar field into a superposition of an electro- magnetic field and dust matter with negative 
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energy density, which turns out to be a methodologically important development and general- 
ization of these models. This allows us to apply methods of Tolman's problem to calculate the 
model (see [20], [21], [22]); these methods were generalized and further developed by Shatskiy 
(see [23] and pi]). 

To further generalize and develop this method, we introduce the cosmological A-term into 
the model. This allows obtaining a principally new solution (see Section [6]) for a static spherical 
model of the Multiverse. This model includes a infinite number of spherical worlds connected 
by throats. To our knowledge, this is the first analytic model of this type. 

In the obtained solution, the Multiverse can have its total energy density positive everywhere 
in space. In addition, this solution can be generalized by the same method to the case of a 
dynamic model shifted from equilibrium by an excess (or shortage) of dust or by the A-term. 

Some methodological details see in [23] 



2 Einstein equations 

We use the Armendariz-Picon static spherically symmetric solution of the Einstein equa- 
tions [1^ for the description of the unperturb WH mode0: 

ds^ = dt" - dR" - r\de^ + sin^ 9 d^^) , r^{R) = + R\ (1) 

The minimum allowed radius of this WH is r^, = g (the radius of the throat^). 

Usually, this solution is related to the energy-momentum tensor of a scalar field (see |19j . 
|17j). In the present paper, we suggest another interpretation of this solution related to a 
different representation of the energy-momentum tensor. This new interpretation of solution 
([T]) can correspond to the energy-momentum tensor represented as 
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The first term in the right-hand side corresponds to the energy-momentum tensor of a static 
magnetic field with the effective magnetic charge q and the second term describes dust matter 
with the negative energy density 
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This possible representation of the tensor is different from its equivalent representation as 
the energy-momentum tensor of a scalar field adopted in previous papers [T7]-[28j. 

An excess (or shortage) of dust (relative to Ed) is not in equilibrium and initiates motion. 
Because dust layers are relatively independent, it is possible to integrate the equations of motion 
for dust in a way similar to the solution of Tolman's problem [20]. In essence, it is the same as 
Tolman's problem in a centrally symmetric and static electric (or magnetic) field for uncharged 
dust. 



^ We use units where c = 1 and G = 1. 

^Here and below, all values at the throat have a subscript 
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We seek the metric tensor of the solution in the form|§: 

ds^ = df - e^dR^ - {de^ + sin^ Odip^), (4) 

where and e'^ are functions of both R and t. Below, we justify this choice of the metric. 

We consider the problem in the presence of a A-term. The Einstein equations corresponding 
to metric (HI) can be written asj: 



SttT/ = 8ne + q^/r^ + A = -e"^ (2rr + r^^ - rr ^A,^) /r^ + (rr ,A,, + r^^) /r^ + , (5) 

SttT^ = gV^^ + A = (2rr,,, + r^) /r^ - e^Vyr^ + l/r^ , (6) 

8vrT,« = = e-^(2r,,-r,AJ/r, (7) 

8nT^ = 8nT^ = -q'/r' + A = - {e"^ [2r,,^ - r ^A,J - 2r „ - rA,, - rAy2 - r,A,} /(2r) . (8) 

Thus, the energy-momentum tensor includes three types of matter: 
a centrally symmetric magnetic field, 
a cosmological A-term, and 
dust matter with the densitjj^ e. 

3 Solution and the Initial Conditions 

Equation ([7]) can be integrated with respect to time: 

e-V=Fi(i?). (9) 
At the throat, the condition r^^ = must be satisfied, and therefore 

i^i(0) = 0. (10) 

By substituting in Eqn Q, we obtain 

1- + Ar^ = (rrl) ^ /r , - Fi + 1 . (11) 
Integrating this equation with respect to time yields 

J - Ar^S + rrl + {1 - Fi)r = F2{R) . (12) 

The functions Fi{R) and F2{R) determine the initial conditions for the velocity and dust density 
distributions. 

With A,( and A^^^ expressed from Eqn ([7]), Eqn (IHl) can be rewritten as 

gV^^ + Ar^ + e'\1^ - r^ - 2rrJ = . (13) 



■^Because dust is pressureless, we can choose the synchronous commoving frame (with the time metric com- 
ponent equal to 1). 

''See, e.g., [31] for the derivation (problem 5 in §100). 

^In the absence of the L term, solution (|4]) coincides with ^ only for the density e equal to Sd [see . ([3])]. 
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When Eqn (E]) is satisfied, integrating Eqn (13) yields an identity. This identity is the 
consequence of setting the time component of metric ([4]) to unity, which justifies 
this choice of the metric. 

It is worth noting from the methodological standpoint that this form of the metric for 
dust is obvious because dust is pressureless and hence the comoving frame is simultaneously 
a synchronous one. But this is not so obvious when a spherically symmetric electromagnetic 
field is added, and the obtained identity is a consequence of the absence of interaction between 
uncharged dust and the electromagnetic field. Moreover, the Lorentz transformation along the 
magnetic field lines does not change this field, which justifies the use of the dust comoving 
reference frame. For motion in this frame (along the field), the electromagnetic field remains 
invariant. If the metric could not be represented in form (jl]), further calculations would be 
impossible. It is this non-obvious fact that allowed us to successfully develop our method and 
to obtain all important results. 

Multiplying Eqn ([5]) by r^r^^ and expressing A^^ from Eqn ([7]), we obtain 

87r£rV^^ + Kr^ - /r -r{l - Fi) - rr"^^^ =0. (14) 
Integrating this with respect to R from to i? with account for ffTOj) yields 

R 

j Sner'^r^^ dR - q^/r + Kr^ /?> - r(l - Fi) - rr^ = -q'^/r^ + Ar^/3 - - r^r'^^^ . (15) 



Here, we took into account that the function r{R = 0) = r^^(t) can be time-dependent. 
Expressing rr^^ from (IT^ and substituting it in (HM . we obtain 

R 

J Snerh,^ dR = F2{R) - gV^o " + Ar,V3 - r^r^^ . (16) 



The derivative of this expression with respect to R yields one more integral of motion (which 
expresses mass conservation in the comoving volume): 

SvrerV, = ^ . (17) 



The Bianchi identified: 



1 d{Tl^^) T"' dg^i 



rpn V fc V J / J"' n 



are contained in the Einstein equations, their components k = t,6,(p vanish identically, and 
the component k = R yields the result in (ITTl) . which is an analogue of the integral for the 
pure-dust solution (see [31], §103) in the theory of the evolution of a dust cloud. The explicit 
form of the function e^{R, t) = r'^^/Fi and an implicit form of the function r^(i?, t) are obtained 
in Section m 

Below, we use the index "i" to denote all quantities at the initial instant t = 0. 

The velocity distribution is set to zero initially: 



^2 



= (19) 

(t=0) ^ ^ 



^Here, g is the determinant of the metric tensor. 
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A nonzero initial velocity is considered in Appendix \M We note that in all cases, the choice of 
initial conditions is not arbitrary and must satisfy constraints imposed by ([5]) and ([7j). These 
equations do not contain second time derivatives and, as is well known (see |27j), under the 
condition that t = const at the initial cross section, they are automatically satisfied for all 
values t by virtue of the other Einstein equations. 

We choose the initial coordinate scale of i? in (jll) such that 



q^ + R\ ed, = edin) = -^. (20) 
where ri{R) coincides with ([1]) and Eqn (fTOj) is automatically satisfied. We set 



R n 
- f 87r{ei - Ed^rln^^ dR = - f Svre.r^ dn + 2^(l- ^) (21) 



The quantity s has the meaning of the doubled Newtonian potential of dust matter excessive 
relative to Ed [see (!20|) ]. The relation of this potential to the functions F\ and F2 follows from 
(fT5ll and (flGll at the initial instant: 

g2 A 



i^i = l-^-^-9(-.^-?Vn), (22) 



M^ + 2g-^). (23) 



Equation ^ implies that Fx > 0. According to this imposes constraints on the distribu- 
tion s{R). 

From Eqn ([6]), we can find r at the initial instant: 



s An Aq^ , 
att = 0: r„ = -- + - + ^ (24) 



Equation flM|) shows that for A = 0, the dynamics are absent in the region were s = (at 
least up to the instant of intersection of the dust layers, when the model considered becomes 
inapphcable). 



4 Solution without the A-term 

We consider the case A = in more detail. 

For A = 0, it follows from Eqn that the WH throat size remains constant \r^{t) = q] 
in the whole region of the allowed values of R. and t for the sought solution, until 
matter starts flowing through the throat. In addition, expression ( 1241) implies that 
the dynamics of matter in a WH depend only on the internal layers of matter (with 
smaller values of r) and are independent of the external layers. Therefore, matter 
on one side of the throat does not influence matter on the opposite side (unless 
there is matter flowing through the throat or layer intersection occurs). 

It can be verified directly from fl24|) that at the initial instant, the acceleration r_jj is negative 
for s > 0. Hence, the excessive mass of dust is to collapse. 
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Figure 1: Diagrams t{ri/q): 1 - for light geodesies; 2 - for r = const; 3 - for the apparent horizon 
and the inner horizon; 4 - for the hmiting stopping time; 5 - for the throat r = q; 6 - for the inter- 
section times of adjacent dust layers; 7 - the dependence s{ri/q) = 0.7(rj/g — 1) ■ exp(2 — rj/g) 
— the scale for s (vertical axis) coincides with the time scale t. 



We introduce the definition of the apparent horizon. 

The criterion of the absence of a apparent horizon during the collapse is given 

by 

V < 1 (25) 

for the dust velocity [see also [23], [29], [30], [21] and f|T2|) ]. where the dust velocity V is determined 

as 

{s + 2q'^/rf)ri/r — s — q^/r'^ — q^/r"^ 
1 — s — q^ /rf 

For g = 0, this result coincides with the solution for collapsing dust (see j231|). 

To avoid a 0/0 ambiguity for the function V"^ at the throat at the initial instant [Fi [r^ r^] 0), 
it is necessary to ensure a more rapid convergence to zero of the function r^^ (t = 0, — r^) than 
the function Fi(rj r^,) (see Appendix |A|) . 



it ' it' 1c. r,2 //r.2 ^ ' 



6 



1.0 1.5 2.0 2.5 3.0 3.5 fi 4.0 




Figure 2: Diagrams t{ri/q): 1 - for light geodesies; 2 - for r = const; 3 - for the ap- 
parent horizon and the inner horizon; 4 - for the hmiting stopping time; 5 - for the 
throat r = q; 6 - for the intersection times of adjacent dust layers; 7 - the dependence 
s{ri/q) = (1 — g/rj)(2 + 0.5(1 — q/ri))q/ri — the scale for s (vertical axis) coincides with the 
time scale t. 

The law of motion of a dust layer can be found from f|T2|) : 




— r dr 



—r dr 



(27) 



t = 



^{Fi + Far - q^ 



I ri\/{l- r/ri)[{s + q'^/rf)r/ri - q^/rf] 



This quadrature can be elementarily integrated: 



t = 



\Jil- r/ri)[{s + q^/rf)r/ri - q^/rf 



s + 2q^/rf 



■arccos 1 — 



2{s + q^/rf){l -r/n) 




2{s + qyr^)y^ 



s 



(28) 
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or in another form 



J (1 - r/ri) \{s + q^/riy/ri - q^/rf] g + 2q^/rf 

' arctg^ 



(s + g2/r?)(l - r/vi) 



s + q^/rf (s + g^/r?)'^/^ ^ {s + q"^ / rf)r / Vi — q"^ / r. 

(29) 

Hence, the motion of dust layers depends on the excessive potential s and the 
initial distance from the WH throat: Vi/q. 

However, this model is physically correct only until a possible intersection of dust layers. 
This is because after the intersection, an incoming energy flux appears in the matter-comoving 
frame (from the intersected matter layers), which has not been taken into account in the energy- 
momentum tensor. The intersection of adjacent matter layers corresponds to an infinite energy 
density. This corresponds to the condition r^^ = (or ^ = 0) [see (l20l) ]. Differentiation of 
f l29l) with respect to allows explicitly determining the function ^(r,s,q). Differentiation of 
the left-hand side of (1291) yields zero. We let a prime denote the derivative with respect to r^, 
express all distances in the units of units r^, and omit intermediate bulky calculations; the final 
expression for this function is given by 



(1 - r) [s{s' - 2g^)(2g^ - r{s + q^}) - q\s + 2q^){2s' + s)] 



[2s{s + g2) + (2g2 _ s'){s + V)] ^(1 - r)[(s + q^)r - q^] 

arctg^ 



(g + g2)(l -r) 

2r(s + q^Yl"^ [s + q'^)r — q^ 



(30) 



It is impossible to use this formula to determine the instant of intersection of nonadjacent layers 
because the intersection with nonadjacent layers does not lead to any irregularities for dust. 
Clearly, the intersection of nonadjacent layers occurs after the intersection of adjacent ones and 
cannot propagate with a superluminal velocity. Therefore, after constructing diagrams t(rj) for 
light cones and the curve corresponding to the intersection of adjacent layers (r' = 0), we can 
determine regions that are definitely free from layer intersection (see the Figure [1] or [2]) . 
Using formula (1261) . we can reach several important conclusions. 

1. The formation of a horizon is possible only for a sufficiently large parameter 

. > 2l f 1 - ^] (31) 



After the appearance of a horizon, the WH becomes "non-traversable". However, 
dust layers can already start intersecting before the horizon is reached. 

2. Values of the function t(r, r^) are bounded by the maximum time tgtop (from 
the beginning of motion until the stop). In (1281) . this time corresponds to the argument 
of arccos being minus unity (or to the zero value of the square root expression in the same 
formula) : 

_ 7rr,(g + 2gVrf) 

2(s + gVr2)3/2 ^-"^^ 

The time tstop corresponds to the deviation Ar ^ srf/q^ from the initial position. For small 
deviations (corresponding to small values of s), harmonic oscillations with the 
period T = 2tstop ~ 27irf/q occur (see Appendix |B|) . 

3. The existence of the second (smaller) root of the equation V'^{r) = 1 implies 
the appearance of the second (inner) horizon in the system (see the Figures [1] and 
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4. For parameter s = 2g(l — q/ri)/ri (or e = 0) the scalar g'''r.ir.k = (1 — q/r)^ (scalar 
equal to the scalar of the extremal Reisner-Nordstrem solution). Therefore: if the energy 
density tends to zero — the solution ( I29h tends to the extremal Reisner-Nordstrem 
solution of black hole (in the comoving, free- falling frame). 



5 Solution for s = with the A-term 

We consider the case where s = 0. Eqn ([6]) with Q and (!22|) then becomes 



(33) 



and Eqn ( 1121) takes the form 



r — r,- 



(r - n) + 



A 



2^2 V ' ' 3 



r + r j + 



rr,- 



(34) 



These two equations imply that at s = 0, matter starts expanding from the rest state (see fTOj) : 
inflation due to the A-term). If the dimensionless parameter 



a = Ag^ 



(35) 



exceeds some critical value acr, the inflation continues unlimitedly in the entire 
volume {q < Ti < r < oo) until an outer horizon forms, V'^{r) = 1. The value is found 
from the condition of the maximum of the expression for the parameter a, which is obtained 
by equating the second factor (in square brackets) in Eqn fl51|) to zero: 



MAX 



3g^(l - r)r'^ 



1 + r + r^g^ 



0.22, r = ri/r, q = q/ri 



(36) 



For a < Ocr, a region appears where the inflation stops. This region starts emerging at 
the throat and extends toward larger values of as the parameter a decreases. 



As r/rj — > oo, the inflation occurs exponentially: r cx exp (ty A/3). 
The function r' can be expressed in quadratures (see Appendix ICj) : 



1 + \/r — Ti exp 



Po(x) dx 



Q2{r) 



P,(r) 



2r2 



2 V r-V, 



a + (r + Ti 



exp / P2{x) dx\ dr , 



rf(r-ri) 



+ a ( r + + ^ 



Q2 



aq^ 
2r^r^ 



rf(r-ri) 



+ a (r + Ti + — 



(37) 



Using this equation, we can find regions where the inflation occurs without dust layer intersec- 
tions. 



6 Model of the Multiverse 

Dust with the positive s initially accelerates toward the center (without the A-term). The 
A-term without excessive dust leads to the original inflationary solution [see (^Mj]- At the WH 
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throat, the potential s = 0, and the A-term at the throat provides nonzero matter acceleration. 
This contribution cannot be compensated at or near the throat [see ([6])], where r^^ = r^^^ = 
and r = q. Hence, there is no static solution for a WH with the throat radius r^^ = q and the 
A-term. 

Nevertheless, a static solution exists for a WH with the A-term and 7^ q. This solution 
can be easily derived from Eqns and (with r^^ = and r = 0). 

The metric of a static WH is determined by expression ([1]), taking into account that 
r\R)^q^ + R\ 

From (6), we obtain r^^ = 1 — g^/r^ — Ar^ and easily find the expression for the throat 
radius r„, 



Ar:-r: + q' = ^ = ^-^^^^^ q\l + Aq') (at Ag^ ^ 0) , (3^ 
and the dependence r{R): 



The distribution of e for this solution is determinecill from 

A g2 



(40) 



Thus, for the static solution with the A-term, the throat radius and dust density are larger 
than without the A-term. 

The total energy density for the static solution with the A-term is 



, 3A g2 

T/ = ^ 41 



The condition for the value T* to be nonnegative everywhere in space (for r > r^) has the form 

T/ > at a = Ag2 > — (42) 

16 

Therefore, in the presence of the A-term, solutions for wormholes with a positive 
total energy density can be found. 

From fl39|) . we can obtain the maximum allowed radius of the static metric with the A-term: 



2 1 -I- yi - 4Ag^ 1 
fmax = ^ ' J (at Ag 0) . (43) 

Beyond this radius, the Universe starts contracting again until a new throat occurs. 

Solution ( 1391) describes a static Multiverse with an infinite number of throats. 
When there is no charge (g = 0), this solution becomes a cosmological solution 
without wormholes, corresponding to a closed isotropic universe (see §112, [31]). 

When there is an excess (or shortage) of the dust part of the energy density or the A- 
term relative to (HUj) . the obtained solution for the Multiverse becomes dynamical. Its analytic 
study is complicated by the need to solve a fourth-order algebraic equation and to calculate a 
quadrature similar to fl37j) for the instant of dust layer intersection. 



^Here, the unknown value r^^^ — /r^ — Ar can be obtained from expression ^ or by directly differentiating 
Eqn 
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7 Conclusion 

In this paper, we have generahzed and developed the method suggested earlier in papers [20], 
[22] and applied it to new problems in modern cosmology. We found and analyzed analytic 
solutions of the general relativity equations describing the dynamics of a traversable wormhole. 
The results obtained are important for the analysis of the general properties of traversable 
wormholes. We also obtained a solution describing a spherically symmetric model of the Mul- 
tiverse. We have not analyzed the properties of geodesies that describe the motion of particles 
and other matter (energy and information) through a wormhole and its vicinity. These prob- 
lems, which are mostly important for the analysis of the possible observational appearance of 
such objects, are considered elsewhere (see, e.g., [TT]). 

To conclude, we remark on the term " Multiverse" . This term is used in physics and cos- 
mology in two different senses. 

First, this term assumes the possibility of the parallel co- existence of many or even an 
infinite number of different worlds, possibly emerged from a quantum vacuum (in some sense, 
in different places at different times). In this paper, this term is used exactly in this sense. 

Second, this term is sometimes used to denote a sample of different realities in the Everett 
interpretation of quantum mechanics. The usage of one term for different notions is sometimes 
confusing. In our opinion, different terms should be used for these concepts. We propose keeping 
the term "Multiverse" for cosmology, and refer to the set of Everett worlds (following the 
proposal by M.B. Mensky) as the "Alterverse", keeping in mind different classical alternatives 
of the Everett world. 
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A Initial conditions with a nonzero velocity 



At the WH throat, the expression for in (!26l) involves a 0/0 ambiguity. To resolve this 
ambiguity, we consider the model with a nonzero initial velocity of matter. We introduce the 
notation 



1 - q/ri = a , ri/r -1=P, 



(i=0) 



= 7 



(44) 



and assume that A = 0. The expressions for the functions Fi and F2 with 7 7^ [see ( fT5|) and 
( IT2l) ] are given by 

(45) 



Fi = 1 - s - gV^j^ + 7 - loQ/ri , F2 = [s + 2gV^»^ + %Q/ri 
With (fT2|) and a nonzero initial velocity, the expression for V"^ in (126!) takes the form 
-s - q^/rf + 7 - 7(,g/ri + (s + ^q^/rf + ^^qlr^rijr - q^/r"^ 
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In the vicinity of the WH throat close to the initial instant, the functions a and /? take small 
values. From fl2Tl) . in the linear order in a, we have 



s ~ , K = 87rr'^(ei — £d) ~ STrSf^q'^ + 2 < 2 . 
We consider Eqn f H6|) in the linear order in a and (3. After all transformations, we finally obtain 



7-7„(l-a) + a(2-fi:) a(7„ + 2 - /t) ' ^ ^ 

It follows that the condition for V^^ to be nonsingular at the throat is that 7^ = 0. Then 

^ J (48) 

7 + a(2-fi:) ^ ^ 

Expression ( HHl) for is regular and has no ambiguities if the function 7 tends to zero at the 
throat faster than the function a. In this case, the rate of change of the radius r vanishes at 
the throat (as must be the case by the definition of the throat). In the linear order in /3, the 
function V^^ is time- independent near the throat. 



B Study of intersections of adjacent layers 

We consider the model without the A-term for s > 0. In this case, the excessive matter 
collapses, i.e., r < r^. 

Lemma. There always exists a nonzero time interval [0,t] inside which no inter- 
sections of adjacent layers occur. 

Equation fl30|) implies that this is obvious for a nonzero potential s. The only point at which 
this statement must be proved is the point corresponding to s = 0. 
We study the asymptotic regime s — >• in more detail: 

s « q^/r1 < 1 (49) 

From the allowed range of the radius (until the stopping point) r > q^ /{sTi + jvi)^ we deduce 

1 - - < —^YTi ^ -TTi « 1 (^0) 
Tj s + q'^/rf q Irf 

Because the last term in fl30l) contains no s in the denominator, it can be neglected compared 
to the previous term. This can be proved in more detail by considering the there exhaustive 
options for the ratio of the numerator and denominator in the square root expression of the 
arctan: 

1. (s + q^ /rf){l — r/ri) « (s + q^ /rf)r/ri — q^ /rf 

2. (s + q^ /rf){l — r/r-^ ~ (s + q^ /rf)r/ri — q^ /rf 

3. (s + q^ /rf){l — r/r-^ » {s + q^ /rf)r/ri — q^ /rf 

Keeping the leading terms, we obtain the asymptotic form of Eqn (130|) : 

, {l-r/n)[-q^-2q^-2ns'] . . , 
r ^ IH ^— ^ = 1 - 2riS (1 - rlri)/s (51) 



12 



October 2, 2008 



Alexander Shatskiy, I.D. Novikov, N.S. Kardashev 



In a similar way, the first term in the right-hand side of expression fl28p can be neglected 
(compared to the second term), and hence asymptotically with respect to time, we have 



tq 



2g2 



cos ^ ^1 2 ~ 



Ti + 



srj 



cos I - 1 



sr-; ■ ■ 2g2 

Comparing Eqns flSTj) and fl52l) . we obtain the sought asymptotic form as s — > 0: 

2rfs' 



(52) 



1 



■ sm 



tq_ 
2r? 



(53) 



This proves the lemma. 

In addition, Eqn ( |52i) implies harmonic dynamics (oscillations) for s — > and 2rfs'/g^ < 1. 



C Quadrature for r' in the solution with the A-term 

We differentiate r^^ with respect to rj: 

dr'^^ , r,dr' dr' 1 rfr^^ 
± = 2r r = 2r =^ = ■ - 

dvj '* '* '* dr dr 2r'^ drj 



(54) 



Substituting Eqn (IMI) in (15^ and differentiating it with respect to rj, we obtain the following 
equation for the function y{r) = r': 

dy 



dr 



+ P{r)y = Q{r) 



(55) 



Here, the functions P(r) = Pi{r) + P2{r) and Q{r) = Qi{r) + Q2{r) are given by (for conve- 
nience, Tj is set to unity below): 

l/(2r2) - (r - l)/r3 - (a/2)(l - q^r^) 



Pi = 

Qi 



2(r- 1) ' 
-1 



P2 = 
, Q2 



(1 — r)/r^ + a(l + r + g^/r) 
l/(2r2) + (r - l)/r2 + (a/2)(l - q^/r) 



2(r - 1) 

Equation (l55l) has a standard solution. 



(1 — r)/r^ + a(l + r + g^/r) 



(56) 



(57) 



The constants Ci, ri and r2 are determined by the initial conditions. The exponential can be 
rewritten as (the first term with Pi is integrated): 





r 






r 




y{r) = exp 


— J P{x) dx 


\ 


J Q{r) ■ exp 


+ J P{x) dx 










^r2 







exp 



P{x) dx 



iri 



Co 



■ exp 



P2{x) dx 



(58) 



After that, the part of integral (157|) corresponding to Qi can also be integrated (by parts), with 
the result 



y{r) = \/r — \ ■ exp 
P2{r) 



P^ix) dx 



exp 



Poix) dx 



1 



r2 



exp 



+ / P^ix) dx 



dr + 



r2 
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We can then redefine the lower integration limit and determine the constant Ci from the initial 
conditions. The final expression for the function y{r) = r' is given in fl37j) . 
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